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Abstract 



We present an unquenched Nf = 2 lattice computation of the Bx parameter which controls 
K° — K° oscillations. A partially quenched setup is employed with two maximally twisted 
dynamical (sea) light Wilson quarks, and valence quarks of both the maximally twisted and 
the Osterwalder-Seiler variety. Suitable combinations of these two kinds of valence quarks 
lead to a lattice definition of the Bx parameter which is both multiplicatively renormalizable 
and 0(a) improved. Employing the non-perturbative RI-MOM scheme, in the continuum 
limit and at the physical value of the pion mass we get B^ Gl = 0.729 ± 0.030, a number 
well in line with the existing quenched and unquenched determinations. 
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1 Introduction 

Strong interaction effects in neutral K° — K° meson oscillations are parametrized by 
the so-called "bag parameter" [1]. Several lattice computations of B^ based on 
different lattice fermion discretizations have been performed in the years, mostly in 
the quenched approximation. Attention has now shifted to unquenched estimates. 
See refs. [2, 3] for an updated review on the subject. A compilation of recent data 
is provided in Fig. 1 where renormalization group invariant (RGI) values of Bk 
(denoted as B^ G1 in this paper) are reported. Remarkably, no significant dependence 
of B^ Gl on the number Nf of sea quark flavours is observed. 
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Figure 1: A compilation of quenched and unquenched results for . Data are 
from refs. [4-12], respectively. The symbol * denotes determinations where the 
continuum extrapolation was carried out. 



The standard way of computing Bk with Wilson fermions yields limited accu- 
racy due to partial control of two sources of systematic effects. 

• First of all, the AS = 2 four-fermion operator relevant for the B^ calculation, 
namely 

AS = 2 = s^(l- 75 )ds lll (l- 75 )d (1.1) 

mixes under renormalization with four other operators of the same dimension, but 
belonging to different chiral representations [13] (the so-called "wrong chirality mix- 
ing" phenomenon) . The reason behind this feature is the lack of chiral symmetry of 
the regularized lattice action. In other regularizations (e.g. staggered or Ginsparg- 
Wilson fermions) partial or full chiral symmetry ensures that the operator (1.1) is 
multiplicatively renormalizable (i.e. does not mix with other operators). An im- 
portant detail here is that mixing actually only affects the parity-even component 
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of the four-fermion operator (1.1). In the following we will denote it by Ovv+AA 
dropping the superscript AS = 2 for simplicity. Its parity-odd component, similarly 
called Ova+AVj is protected against mixing by discrete symmetries, irrespective of 
the status of chiral symmetry in the lattice action, as shown in [14, 15]. 

• Secondly while a lattice B>k estimate, based on staggered or Ginsparg- Wilson 
quarks, would be affected by only 0(a 2 ) cutoff artifacts at finite gauge coupling, 
computations based on plain Wilson fermions suffer from 0(a) discretization errors. 
A well known remedy to this state of affairs is to make use of the Symanzik im- 
provement strategy which entails the inclusion of the Clover term in the action, as 
well as a number of dimension-7 counter-terms in the operator. However the non- 
perturbative determination of the coefficients with which the counter-term operators 
should enter would greatly complicate the computation and increase the systematic 
error. 

A proposal to circumvent the first of these problems was put forward in ref. [16] 
where by the use of Ward-Takahashi identities (WTIs) it was shown that the ma- 
trix element {K°\O V y +A a(0)\K ) can be expressed in terms of the matrix element 
f d 4 x(K°\OyA+Ay(0)P(x)\K }, where P is the pseudoscalar quark density. It is seen 
that the latter involves the multiplicatively renormalizable parity-odd part of the 
four-fermion operator (1.1). The price to pay for this result is that now one has to 
compute a four-point correlation function (while normally B>k is obtained from the 
large-time asymptotic limit of a three-point correlation function). This implies in- 
creased statistical fluctuations which in practice may offset the gain stemming from 
the absence of wrong chirality mixings [17]. 

Other attempts to avoid this problem have been tried out in refs. [5,18], inspired 
by the twisted-mass formulation of lattice QCD (tm-LQCD) [19]. Two possibilities 
have been explored. 

The first one consists in using a lattice fermion action with a maximally twisted 
up-down fermion doublet and a standard (untwisted) Wilson strange quark. As a 
second possibility, viable only in the quenched approximation (for lack of reality of 
the resulting fermion determinant), a degenerate down-strange doublet with twist 
angle 7r/4 is introduced. In both variants the change of variables that brings the tm- 
LQCD quark action to the standard Wilson form happens to map the parity-even 
AS = 2 operator onto its parity-odd counterpart, thereby implying multiplicative 
renormalizability. 

All the methods described above [5, 16, 18], while avoiding the operator mix- 
ing problem, lead to B>k estimates which suffer from 0(a) discretization effects. 
Achieving 0(a) improvement, without the use of Symanzik counter-terms, necessi- 
tates having all fermions at maximal twist (Mtm-LQCD [20]). Clearly, this must be 
done while preserving the mapping from parity-even to parity-odd AS = 2 operator, 
so as not to loose multiplicative renormalizability (see the discussion in ref. [21]). 

This can actually be achieved by using a different regularization for sea and 
valence quarks. The idea proposed in ref. [22] and adopted in the present study is 
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to have maximally twisted sea quarks in combination with Osterwalder-Seiler (OS) 
valence quarks. We recall that any OS fermion can be thought of as a component 
of a maximally twisted doublet with twist angle +7r/2 or —ir/2, depending on the 
value of its Wilson r-parameter. 

Then the relevant AS = 2 four-fermion operator is defined in terms of four 
distinct (maximally twisted) valence flavours, three with twist angle, say, +vr/2, and 
the fourth with twist angle —tt/2. The discrete symmetries of the valence quark 
action combined with the structure of the resulting four-fermion operator ensure 
that the latter is multiplicatively renormalizable, while its K° — K° matrix element 
remains automatically 0(a) improved [22]. The same is then true for Bk- 

In the setup of ref. [22] unitarity violations occur as a consequence of the differ- 
ent regularization of sea and valence quarks. However, if renormalized masses of sea 
and valence fermions are matched, one can prove [22] that unitarity violations are 
mere 0(a 2 ) effects. Moreover, it just happens that the valence quark content of the 
ivT°-meson consists of a strange/down pair with, say, the same twist angle, while the 
^T°-meson has necessarily a quark-anti-quark pair with opposite twist angles. Thus 
the two pseudoscalar mesons have masses that at non-zero lattice spacing differ by 
(numerically important) 0(a 2 ) effects, which mainly come from the lattice artifacts 
(non-vanishing in the chiral limit) of the i^°-meson mass. The ET°-meson mass on 
the other hand exhibits only 0(a 2 (/i^ + /x s )) discretization errors (here fj,£/ s denotes 
the mass of the light /strange quark). This effect has been recently studied in the 
quenched approximation adding the clover term in ref. [7], where it is numerically 
demonstrated that the scaling behaviour of both Bk and the decay constants of the 
K° and K° mesons is only weakly affected (see also below). The situation about 
the 0(a 2 ) mass-splitting between K° and K° is similar to the one already met in 
the unitary Mtm-LQCD framework, where a large 0(a 2 ) artefact is only observed 
in the mass of the neutral pseudoscalar meson as it is not associated to a conserved 
lattice axial current [23,24]. 

In this paper we present an unquenched computation of Bk based on the strat- 
egy of ref. [22] which exploits the lattice data coming from the state-of-the-art tm- 
LQCD simulations carried out by the ETM Collaboration (ETMC) with N f = 2 
dynamical (sea quark) flavours, at three lattice spacings and "light" pseudoscalar 
meson masses in the range 280 MeV < mps < 550 MeV. Renormalization is carried 
out non-perturbatively in the RI-MOM scheme [25]. In the continuum limit and at 
the physical value of the pion mass we get 

Bf^ 1 = 0.729(25)(17) e> 0.729(30) @ N f = 2 , (1.2) 

where in the first expression the two errors quoted are the statistical one (0.025) and 
(our estimate of) the systematic uncertainty (0.017, resulting from the quadratic sum 
of 0.004 from renormalization, 0.009 from control/removal of cutoff effects and 0.014 
from extrapolation/interpolation to the physical quark mass point). In the second 



3 



expression the total error is obtained by a sum in quadrature of the statistical and 
systematic ones. The error budget is further discussed in sect. 3.2. The specification 
@ Nj = 2 is to remind that the whole computation of B^ Gl has been carried out in 
QCD with two light (u and d) dynamical quark flavours. The quality of this result 
is extremely satisfactory and fully comparable with other quenched and unquenched 
numbers (see Fig. 1 for a compilation of recent determinations). 

The outline of the paper is as follows. In sect. 2 we briefly recall the theoretical 
framework which allows to obtain an 0(a) improved and multiplicatively renormal- 
izable lattice expression of the K° — K° matrix element of the four-fermion AS* = 2 
operator (1.1). In sect. 3 we give the details of our numerical analysis. Conclusions 
can be found in sect. 4. We defer to a couple of Appendices a few more technical 
considerations. Preliminary reports about the present work have already appeared 
in refs. [3,26,27] 



2 Twisted and Osterwalder— Seiler valence quarks at maximal twist 

The lattice setup of our study is that of ref . [22] . The regularization of sea quarks is 
different from that of valence quarks. The former is a standard tm-LQCD regulariza- 
tion with a doublet of degenerate, maximally twisted Wilson fermions (representing 
up and down flavours). In the so-called "physical" basis, the Mtm-LQCD fermionic 
action takes the form 

<S M tm = a 4 X^( V M + V*) - i 7 5T 3 [M cr - | J] V M V ^] + Vsea}?P(x) , 

(2.1) 

where the Wilson's r parameter has been set to unity, tp(x) is the quark flavour 
doublet, and are nearest-neighbour forward and backward lattice covariant 
derivatives, fi sea is the (twisted) sea quark mass and M cr the critical mass. The 
lattice pure gauge action is the tree- level improved action of ref. [28], routinely 
adopted by the ETM Collaboration [29] in studies with two light sea quark flavours. 
Valence quarks are regularized as Osterwalder-Seiler (OS) fermions [30]. The action 
of each OS valence flavour qj reads 

Sos = a 4 X ?/(*){ \ E 7^(V M +V*)-i75r/ [M cr -f £ V* V„] + N )q f (x) , (2.2) 

x n fi 

where aM cr is the same number as in eq. (2.1). The critical mass M CI is non- 
perturbatively tuned to its optimal value [31] as described in ref. [29]. This ensures 
0(a)-improvement of physical observables and control of the leading chirally en- 
hanced cutoff effects 1 . 

x At the formal level a ghost action is also involved to cancel the fermionic determinant, arising 
from the integration over the qf and q~f valence degrees of freedom [22]. The presence of ghosts 
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2.1 Lattice operators 

For the computation of interest, it is convenient to introduce four OS valence quark 
flavours q±, q 2 , Q3, </4- Eventually q± and q$ will be identified with the strange quark 
by setting fi± = fi^ = fih, and qi and q^ with the down quark by setting fi2 = 1^4 = He- 
The key point is that the four-fermion operator 

Qw+AA = 2 {[ft 7^2] fe7^4] + [<?17m75<?2] [937^7594] + <-> ft)} (2.3) 

is multiplicatively renormalizable once the Wilson parameters appearing in the ac- 
tion of the OS fermions are taken to satisfy the relations 

n = r 2 = r 3 = -r 4 . (2.4) 

The factor 2 in the r.h.s. of eq. (2.3) guarantees that the correlator (2.11), where 
Qw+AA is inserted between the interpolating K° = q 2 q\ and K° = q^q^ fields, is 
normalized as in continuum QCD (as one can check by direct application of the 
Wick theorem). The proof of the multiplicative renormalizability of the operator 
Qw+AA is based on the discrete symmetries of the OS action and is provided in 
full detail in ref. [22]. A simpler proof is readily obtained in the so-called "twisted 
basis". In fact, upon performing the chiral field rotation 

q f -+ qf = e^/75/4 q ^ 

q f -+ qf = q je ™n*l* (2.5) 

from the "physical" (qf) to the "twisted" (qf) basis, the valence quark action in 
the massless limit (all fi/s set to zero) takes the standard (untwisted) Wilson form, 
while the parity-even operator Qw+AA gets rotated onto the parity-odd Qva+av- 
The latter, being expressed in terms of quark fields with standard Wilson action, 
is known [14, 15] to be multiplicatively renormalizable. Moreover, the argument 
implies that the renormalized operator is given by the formula 

[Qw+aaJr = Zva+av Qw+aa , (2.6) 

where, for consistency with literature on Wilson fermions, the renormalization con- 
stant is named after the form (VA + AV) the operator takes in the "twisted" basis. 
In actual simulations we fixed the common value in eq. (2.4) by setting n = 1. 
In order to compute Bk, we also need to consider the axial currents 

A ^ = Qil^5Q2 A 3 * = <?37m75«4 (2.7) 

which, as a consequence of the relations between the Wilson parameters r\, . . . 
specified in eq. (2.4), when passing from the "physical" to the "twisted" basis take 



has no practical consequences in the evaluation of correlation functions of operators where no ghost 
fields appear. 
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the form of axial and vector operators, respectively. Consequently they renormalize 
according to the relations 



[Af] K =Z A Af [Af] K =Z v Af. (2.8) 

2.2 Lattice correlation functions 

Next we proceed to construct the correlation functions, involving the operators that 
are necessary to evaluate Bk- The lattice is taken to be of size L 3 • T, with periodic 
boundary conditions on all fields and in all directions, except for quark fields in the 
time direction which satisfy antiperiodic boundary conditions. At a reference time 
slice yo, we define a 11 K- meson wall" with pseudoscalar quantum numbers and q~2 
and qi quark fields, namely 

W 21 (y ) = (|) 3 J>2(y,yo)759i(y,yo). (2.9) 
y 

A second "i^-meson wall", W A3 (yo + T/2), with q^ and q^ quark fields, is placed at 
the time slice yo + T/2. The four-fermion operator, Qvv+AA, is inserted at position 
x = (x,xo), with xo in the range yo x o *C yo + T/2. In this way lattice estimators 
of the bare l?x-parameter can be calculated at several values of xo from the ratio 

involving the correlation functions 2 

C 3 (x ) = ^Y^2(W 4;i (y + ^)Q vv+AA (x)W 2 \yo)), (2.11) 

C2(x ) = (f) 3 ^(A 2 (x)^ 21 (y )}, (2.12) 

X 

C' 2 (x )=Q 3 J2(W^(yo + ^AlHx)). (2.13) 



The signal-to-noise ratio is greatly enhanced by summing over the spatial posi- 
tion of the quark and antiquark fields in each "K-meson wall" as well as summing 
over that of the four-fermion operator. The xo-behaviour of the i?x-estimator ex- 
cludes significant contaminations from excited states (see also sect. 3.1). On the 
other hand in order to reduce the cost of computing the quark propagators we chose 



2 In order to simplify the notation, we do not display the dependence on the reference time slice 
yo- The latter is chosen randomly on each gauge configuration so as to reduce the magnitude of the 
autocorrelation time of Bk lattice estimators. 
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to employ a stochastic method. This consists in computing quark propagators by 
inverting the lattice Dirac operator for the relevant OS valence quark flavours on 
random sources of the form 

S ai S X0 , Z0 S ik S s ^ Z0 \k,z), 7 = 1,2,3,4, z = y ,y + T/2, (2.14) 

with free indices a (spin) , xo (time) , i (colour) , x (space) , while 7 and zq are fixed as 
indicated. The ^-valued 77 vectors carry only colour and three-space indices (colour 
and three-space "dilution" ) and are normalized according to 

( v ^(k,zW^(k',z')) = 5 kk ,5^, . (2.15) 

Given the pattern of Wilson r-parameters specified in eq. (2.4) and the invariance 
under 75-hermitian conjugation combined with rj — V —rf of the OS lattice Dirac 
operator, it is enough to compute for each gauge configuration the propagator of the 
fields q\ and q<i on the random source at yo, an d the propagator of the fields q% and 
^4 on the random source at yo + T/2. The conditions (2.15) guarantee that unbiased 
estimators of the correlators C3, C2 and C 2 are obtained 3 . We found that, for an 
ensemble of a few hundred gauge configurations, employing one set of ^-valued 
random sources (namely those in eq. (2.14)) per gauge configuration is sufficient to 
achieve a good statistical precision for our S^-estimators. 



2.3 Lattice estimator of (renormalized) Bk 

As we said, the masses of (up/down) sea quarks and the masses [^2]r, [i^aIr of 
(down) valence quarks must be tuned to the same (physical) value. At the same 
time the valence mass parameters [/j,i]r and [/j,s]r should be adjusted to the strange 
quark mass in order for the external states to be identified with if-particles 4 . In 

this way the renormalized ratio R{xq) = — -R(^o) computed with distinct OS 

valence flavours in the partially quenched set up specified above, will differ from its 
continuum limit by only 0(a 2 ) discretization errors. 

One peculiar feature of our approach which we have to keep in mind in the anal- 
ysis of simulation data is the fact that kaon and anti-kaon masses are not degenerate 
as the two mesons involve differently regularized valence quarks. In our setting the 
kaon is made up of the valence quarks q2 — > qd and q\ — > q s with r\ = r 2 , while the 
anti-kaon consists of the valence quarks q^ —> qd and q% — > q s with r% = —r^. The 
K° - K° mass difference M 12 - M 34 M K - M R is an 0(a 2 ) effect which has an 
origin similar to that of the well known mass splitting one encounters in the pion 

3 This setup was first presented in [26]. 

4 We recall that, since the quark mass renormalization constant of OS valence quarks is indepen- 
dent of the sign of the corresponding Wilson r-parameter [22], at the level of bare masses we must 
simply require n — ^2 = and /ii = [13. 



7 



sector [24]. Naturally the relative effect is much less important in the kaon case as 
the kaon mass is substantially larger than the pion mass 5 . 

Taking this mass difference into account, we write for the three-point correlation 
function the large time expansion (yo *C xq *C yo + T/2) 

L 6 C 3 (x ) -> (0|^ 43 (0)|P 34 )(P 34 |Q yv+ AA(0)|P 21 )(P 21 |^ 21 (0)|0) 

X AM^m ex P[~ Ml2 ( x o - Vo)] exp[-M 34 (y + | ~ so)] ■ (2-16) 

In this expression |P lJ ) denotes a zero three- momentum pseudoscalar state with 
quark content % and (i, j = 1, • • • ,4) and mass M' 3 . Similarly we have 

L 3 C 2 (x ) -+ ( |<(0)|P 21 )(P 21 |^ 21 (0)|0)^exp[-M 12 (xo- yo )], (2.17) 

L 3 C 2 (* ) (0|^ 43 (0)|P 34 ) (P 34 |A 34 (0)|0) ^_ exp[-M 34 (y + ? - *„)] . (2.18) 

Furthermore, defining the pseudoscalar decay constants F 13 through the equations 

(0|^ 2 (0)|P 21 ) = M 12 F 12 (P U \A 3 4 (0)\0) = M 34 P 34 , (2.19) 

one finds that they are not equal, but again they differ by 0(a 2 ) effects. 

The renormalized lattice Bk parameter which we will have to extrapolate to 
the continuum limit is finally obtained by averaging over the xo-plateau (which lies 
between yo and yo + T/2) of the estimator R(xq), viz. 

8 R ^R(r)- ZVA+AV P(~ ) - ZVA+AV (P 3A \QVV + AA(0)\P 21 ) , 99m 

-B KM <= R(x ) = R(x ) - Ml2jpl2M 3 4F 34 > ( 2 - 20 ) 

where the exponentials in eqs. (2.16), (2.17) and (2.18) exactly cancel out and the 
numerical factor 8/3 in the l.h.s. of this equation is there for conventional reasons, 
i.e. it ensures that Br represents the four-fermion matrix element in units of the 
vacuum saturation approximation value. 



2.4 Scaling checks 

In this section we present some tests aimed at checking the size of the discretization 
errors affecting the quantities that enter our computation of Bk (i.e. those that 
appear in eq. (2.20)). In order to keep the issue of lattice artifacts separated from 
those related to the extrapolation to "physical quark mass point" , the present scaling 
tests have been performed at fixed reference values of the renormalized "light" and 
"heavy" (would be strange) quark masses in the MS-scheme at the scale of 2 GeV. 
For convenience we take \x) ~ 40 MeV and jl* h ~ 90 MeV. 

5 The scaling of the K° — K° mass splitting has been studied in detail in the quenched approxi- 
mation in ref. [7]. 
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Figure 2: Data for (M 34 // ) 2 (panel (a)), F 34 //o (panel (b)), A M (panel (c)) and 
Ap (panel (d)) as functions of (a/o) 2 at the reference quark masses /t| ~ 40 MeV 
and fit ~ 90 MeV (in the MS scheme at 2 GeV). We display the best linear fit in 
(a/o) 2 and the corresponding continuum limit result. 



In Fig. 2, panels (a) and (b), we show the scaling behaviour of (M 34 //o) 2 
and F 34 //o, i.e. the mass squared and the decay constant of the "kaon" made up 
of valence quarks and 03, normalized by the chiral limit pion decay constant 
fo ~ 121 MeV determined in ref. [35]. Owing to the choice r% = —r^ the quantities 
M 34 and / 34 are our lattice estimators of the "kaon" mass and decay constant that 
should exhibit best scaling properties. The plots in panels (a) and (b) fully confirm 
this expectation, showing a nice a 2 -scaling with cutoff effects for the lattice "kaon" 
mass and decay constant estimates that are at the level of few percents (from the 
largest lattice spacing to the continuum limit). 

In Fig. 2, panels (c) and (d), we show instead the scaling behaviour of the 
difference of the "kaon" mass squared and decay constant lattice estimators entering 
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our Bk computation, namely 



A M = 



(M 12 ) 2 - (M 34 ) 2 



A F = - 



pl2 _ 



(2.21) 



i?34 



The quantities (2.21) are by construction mere 0(a 2 ) lattice artifacts. We see from 
the panels (c) and (d) that both Am and Ap extrapolate nicely to zero in the 
continuum limit, as expected. The large values of Am in Fig. 2(c) signal the presence 
of a substantial cutoff effect on M 12 , with M 12 on the coarsest lattice differing by 
about 30% from its continuum limit value (~ 570 MeV). On the other hand, the 
lattice artifact difference between F 12 at the largest lattice spacing and its continuum 
limit value is only about 5%. This is readily inferred by comparing Fig. 2(b) and 



Having taken r\ = r-2, the large and positive lattice artifact observed in M 
was expected since for OS valence quark doublets, (unlike the case of twisted-mass 
quark pairs) no isospin component of the lattice isotriplet axial current is conserved 
in the limit of vanishing quark mass. Relying on arguments similar to those given 
in ref. [24], one can however conjecture that the large discretization error in M 12 
is due to dynamically large matrix elements entering the Symanzik description of 
the lattice OS pseudoscalar meson mass, rather than to large coefficients in front of 
some terms of the Symanzik local effective action for OS valence fermions 6 . This 
conjecture suggests that the large cutoff effect detected in the lattice OS pseudoscalar 
meson mass is peculiar to this quantity (or to others directly related to it) but not 
a general feature of physical observables built in terms of OS valence quarks. The 
small cutoff effects observed in F 12 is well in line with such an expectation. 

In Fig. 3 we show the scaling behaviour of jlV) evaluated according to 

eq. (2.20) using two different methods (referred to as Ml and M2) for the evaluation 
of the renormalization constants Zy^+^y and Za- The two methods (discussed in 
sect. 3.2.1 and Appendix A) differ only in the way one deals with 0(a 2 ) artifacts and 
other unwanted systematic effects pertaining to the RI-MOM scheme computation 
of renormalization constants. The somewhat different slope in a 2 of the data-points 
in Fig. 3 should be ascribed to these effects. Nevertheless in both cases the a 2 -scaling 
behaviour of S^^[ t (/i|, /}£) is very good and the extrapolated continuum limit values 
agree very well. Cutoff artifacts in B#fil ct (p>i, A/1) are about 10% with the method 
Ml and a bit larger if the method M2 is adopted. 

Two remarks are in order here. First of all, our choice of the lattice normal- 
ization factor (M 12 F 12 M U F U )~ 1 in the r.h.s. of eq. (2.20) turns out to be very 

6 When extending the analysis of ref. [24] to the OS doublet case one finds that, in contrast to 
the situation in the twisted mass doublet case, matrix elements with both one insertion of the d — 6 
Symanzik local effective action term and two insertions of the d — 5 term contribute to the cutoff 
effects in (M 12 ) 2 . 



Fig. 2(d). 
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Fi gure 3: -B^f^t as function of (afo) 2 , at the same reference quark masses, /2| and 
jit, as in Fig. 2, is shown together with the best fit linear in a 2 and the corresponding 
continuum limit result. The two sets of data, labelled as Ml and M2, come from 
different procedures for evaluating the renormalization constants Z^y+VA and Za, 
expected to agree in the continuum limit. 




beneficial as it partially compensates the discretization errors coming from the ma- 
trix element {P u \Q vv+A a(0)\P 21 }- Secondly, the cutoff effects in the RI-MOM 
scheme computation of Za and Zva+av are significantly reduced when 1-loop per- 
turbative lattice artifacts (i.e. 0(a 2 (7 2 ) effects) are subtracted. The latter have been 
computed in refs. [32,33]. 

In conclusion, in spite of the occurrence of a substantial 0(a 2 ) artifact in the 
mass of the lattice state \P 21 } and the ensuing fact that the {P 34 \Qvv+aa(0)\P 21 ) 
matrix element is evaluated with an 0(a 2 ) four-momentum transfer, the scaling 
of £#iat(Al> Aft.) with a 2 is found to be well under control, thereby allowing for a 
reliable continuum limit extrapolation. The same is thus expected, and found to be 
true (see sect. 3.2.2), also upon approaching the physical quark mass point. 



3 Simulations, data analysis and results 

The ETM Collaboration has generated Nf = 2 configuration ensembles at three 
values of the inverse bare gauge coupling (/3) and at a number of light quark masses 
(Usea)- Several quantities entering the data analysis of this paper, such as the charged 
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pion mass aM^f 1 7 , the low-energy constants afo and clBq and the chiral limit value 
of a/ro, are derived from ETMC data [29,34,35]. In this paper the symbol Q means 
that the quantity Q is renormalized in the MS scheme at the 2 GeV scale. We note 
that, since in our analysis the continuum limit value Bq = ZpBq is employed, in 
order to evaluate the RGI quantity B^i = Bofte one must know fi£ = Zp 1 ^ at 
the /3-values of interest. The necessary renormalization constants of quark bilinears, 
such as Zp, Za and Zy, have been computed non-perturbatively in the RI-MOM 
scheme [25] as reported in ref. [32] and converted to the MS wherever necessary. 
The calculation of the four-fermion operator renormalization constant is new and 
discussed below in subsect. 3.2.1. 

3.1 Extracting bare Bk estimates from lattice data 

In the present computation we deal with hi pseudoscalar mesons ("kaons") made 
of two valence quarks, and q£, with bare masses (/x^, fJ-i)- The values of the mass 
parameters (a/ih, aw) are chosen so as to have 11 pseudoscalar mesons (charged 
"pions" with r u = —r^) with mass in the range 280-550 MeV and hi pseudoscalar 
mesons with mass in the range 450-700 MeV, depending also (see sect. 2) on whether 
the lattice mass M34 <-)■ M^-, or rather Myi Mk, is considered. The value of the 
light quark mass parameter, afi£, is common for sea and valence light (the would-be 
u/d) quark flavours, while the heavier quark (the would-be s) is quenched. With 
this choice of mass parameters we will eventually be able, as discussed in sect. 3.2, 
to make contact with the physical K° and K° mesons by interpolating our results 
in \ih — > n s and extrapolating them in the m — > fi u / d and continuum limit. 

Simulation details are gathered in Tables 1, 2, and 3, where we report, at each 
gauge coupling and for all the available (/i^, //^-combinations, the values (in lattice 
units) of the measured pseudoscalar meson masses and decay constants as well as 
the bare .B^-parameter. The latter (up to a trivial factor 8/3, see eq. (2.20)) is 
denoted by R and it is obtained by averaging over a plateau in r = xo — yo of the 
lattice estimator R(xq) given in eq. (2.10). In the table captions we specify the 
plateaux intervals [(xq — yo) mm /a, (xq — yo) ma,x /a\. Indeed, the quantity R{xq) must 
be evaluated at large time separations. In our case this means requiring that xq (the 
time coordinate where the four-fermion operator is located) is sufficiently far away 
from yo + T/2 and yo (with yo, yo + T/2 the time coordinates of the two "K-meson 
walls" , see sect. 2.2). In this way the K°- and the -K^-state dominate the three-point 
correlator C^{xq). Actually, due to the finite T extension of the lattice, it is also 
necessary that the contribution to the quantum-mechanical representation of C^{xq) 
from states wrapping around-the- world, such as that from |7T7t) or the a 2 -suppressed 
one from \tt°), be negligible compared to the contribution from the vacuum state. 

7 We use this notation to make clear that this is the mass of the "charged" pion which is known 
to be affected by only small 0(a 2 m, a 4 ) lattice artifacts [23,24,31]. 
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/3 = 3.80, a ~ 0.10 fm 
L 3 x T = 24 3 x 48 a 4 


ape an h aM' M aM vl aF' M aF V2 x -plat. 
0.0165 0.2558(08) 0.3393(25) 0.0894(04) 0.0883(15) 
0.0080 0.0200 0.2731(07) 0.3532(23) 0.0913(04) 0.0895(15) [11, 23] 
0.0250 0.2961(07) 0.3718(21) 0.0936(04) 0.0909(15) 


R xo-plat. 
0.627(4) 

0.636(3) [11, 14] 
0.647(3) 


-^meas 

170 


0.0165 0.2712(04) 0.3508(16) 0.0924(03) 0.0900(16) 
0.0110 0.0200 0.2877(04) 0.3644(14) 0.0942(03) 0.0910(16) [11, 23] 
0.0250 0.3098(04) 0.3828(12) 0.0966(03) 0.0924(16) 


0.632(4) 

0.640(4) [11, 14] 
0.651(3) 


180 



Table 1: Pseudoscalar masses, decay constants and the (bare) ratio R at /3 = 3.80. 
The time-plateaux "xo-plat." in column 7 correspond to the interval a -1 fx™ 111 , x™^] 
over which the correlators C2 and C' 2 have been taken, and in column 9 to the 
interval a _1 [(xo — Uo) mm , (xo — yo) max ] over which the correlator C3 has been taken 
(see sect. 2.2). iV meas denotes the number of (well separated) gauge configurations 
on which the various correlators have been evaluated. 



For this purpose it is important that T is large enough with respect to the inverse 
lattice pion mass(es). 

As one can appreciate from Fig. 4, panels (a), (b) and (c), the quality of our 
signal for R(xq) is quite good. Wider plateaux, as well as stronger suppression of 
unwanted vacuum-sector states wrapping around-the- world, are obtained when the 
time extension of the lattice is increased by a factor 4/3 going from a 24 3 x 48 to a 
32 3 x 64 lattice at f3 = 3.9 (i.e. up to T ~ 5.8 fm). This is illustrated in Fig. 4 panel 
(d) where the values of R(xq) and its plateau for the larger lattice is displayed. 
As shown in Fig. 5, no systematic effect in the plateau value of R is visible at 
(3 = 3.90 upon comparing results from the lattices 24 3 x 48 and 32 3 x 64. Similarly 
by comparing data for C^xo) {C' 2 (xq)) from the two lattices, we could also establish 
(see Table 2) the absence of significant finite volume effects in the quantities M12 
and F12 (M34 and F34). In particular we checked that the marginally significant 
finite size effect that one can notice in 0M 34 at /3 = 3.9 and am = 0.0040 has a 
negligible impact on our determination of B^ Gl . 

3.2 Computing B^ G1 at the physical point 

In order to arrive at the RGI value of Bk at the physical point we need to go through 
the following steps: operator renormalization at each lattice spacing, continuum 
(a — > 0) and chiral extrapolation (jii — >■ fi u /d) with fih set to the strange quark mass 
fi s as determined by using the experimental value of the kaon mass. 
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(3 = 3.90, a ~ 0.09 fm 
L 3 x T = 24 3 x 48 a 4 


am an h aM' M aM vl aF' M aF V2 x -plat. 

U.UluU U.iUUUI UJ 1 U.iUOyi 11 1 U.UIi'llUOl U.U IUOIUjI 

0.0040 0.0220 0.2401(05) 0.2915(11) 0.0757(03) 0.0725(09) [11, 23] 
0.0270 0.2619(05) 0.3096(11) 0.0777(03) 0.0759(09) 


R xo-plat. 
n ^s^f^ 

U.uOul O J 

0.608(4) [11, 14] 
0.621(4) 


400 


0.0150 0.2179(08) 0.2762(16) 0.0755(05) 0.0736(10) 
0.0064 0.0220 0.2506(07) 0.3028(15) 0.0785(05) 0.0759(09) [11, 23] 
0.0270 0.2717(07) 0.3204(14) 0.0805(04) 0.0774(09) 


0.602(7) 

0.620(6) [11, 14] 
0.631(6) 


200 


0.0150 0.2283(07) 0.2849(17) 0.0773(03) 0.0755(09) 
0.0085 0.0220 0.2598(07) 0.3109(15) 0.0804(03) 0.0779(08) [11, 23] 
0.0270 0.2803(07) 0.3281(15) 0.0823(03) 0.0794(09) 


0.606(6) 

0.626(5) [11, 14] 
0.637(5) 


200 


0.0150 0.2351(07) 0.2892(14) 0.0784(04) 0.0761(09) 
0.0100 0.0220 0.2659(07) 0.3154(12) 0.0815(04) 0.0787(08) [11, 23] 
0.0270 0.2860(06) 0.3328(12) 0.0834(04) 0.0802(09) 


0.606(8) 

0.625(8) [11, 14] 
0.637(7) 


152 


L 3 xT = 32 3 x 64 a 4 


0.0150 0.2041(04) 0.2644(15) 0.0727(03) 0.0702(11) 
0.0040 0.0220 0.2381(04) 0.2917(17) 0.0758(03) 0.0722(10) [11, 31] 
0.0270 0.2599(04) 0.3096(15) 0.0777(03) 0.0753(09) 


0.592(5) 

0.615(5) [11, 22] 
0.630(5) 


160 


0.0150 0.1982(04) 0.2558(13) 0.0720(03) 0.0701(07) 
0.0030 0.0220 0.2329(04) 0.2838(11) 0.0750(03) 0.0720(08) [11, 31] 
0.0270 0.2550(04) 0.3021(11) 0.0770(03) 0.0745(08) 


0.585(4) 

0.606(4) [11, 22] 
0.619(4) 


300 



Table 2: Same as in Table 1 for f3 = 3.90. 



P = 4.05, a ~ 0.07 fm 
L 3 x T = 32 3 x 64 a 4 


am a\i h aM' M aM vz aF' M aF vz x -plat. 
0.0120 0.1602(08) 0.1931(18) 0.0564(03) 0.0558(07) 
0.0030 0.0150 0.1751(08) 0.2053(17) 0.0578(03) 0.0566(07) [14, 31] 
0.0180 0.1889(08) 0.2169(16) 0.0591(03) 0.0573(07) 


R xo-plat. 
0.560(6) 

0.575(6) [14, 19] 
0.588(6) 


i^meas 

190 


0.0120 0.1739(06) 0.2034(11) 0.0600(04) 0.0585(08) 
0.0060 0.0150 0.1877(06) 0.2153(11) 0.0613(03) 0.0596(08) [14, 31] 
0.0180 0.2007(06) 0.2266(11) 0.0625(03) 0.0605(08) 


0.584(7) 

0.595(7) [14, 19] 
0.605(6) 


150 


0.0120 0.1840(05) 0.2127(09) 0.0615(04) 0.0604(12) 
0.0080 0.0150 0.1972(05) 0.2242(09) 0.0627(04) 0.0616(12) [14, 31] 
0.0180 0.2097(05) 0.2351(09) 0.0638(04) 0.0626(12) 


0.600(6) 

0.609(5) [14, 19] 
0.618(5) 


220 



Table 3: Same as in Table 1 for f3 = 4.05. 
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,8 = 3.80 



(n/if.afu) = (0.0080,0.0200) 
(aw,HA) = (0.0110,0.0250) 



n; 0.70 



''ham 
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Figure 4: Data and time-plateaux for R(xo) plotted vs 2t/T = 2(xq — yo)/T for two 
different combinations of light and heavy quark masses: (a) results at (3 = 3.80; (b) 
results at /3 = 3.90 for the lattice 24 3 x 48; (c) results at f} = 4.05; (d) results at 
f3 = 3.90 for the lattice 32 3 x 64. 



3.2.1 Computing Zav+VA 

Renormalization of Bk requires in particular the calculation of the four-fermion op- 
erator renormalization constant, Z^v+VAi m the mass independent scheme. This 
task was carried out using the RI-MOM approach. We followed the general strat- 
egy presented in ref. [32] employing the two procedures referred there as Ml (or 
extrapolation method) and M2 (or p 2 -window method). In this way one obtains for 
each f3 two estimates of the renormalization constant at the lattice reference scale, 
which we will denote as Z^'^y A (a~ 2 ) , j = 1,2. Using the known NLO anomalous 
dimension of the four fermion operator in the RI-MOM scheme, the RGI quantities 
-^AV+vvi are evaluated. We note that the cutoff effects on the lattice estimators of 
Zav+VA are significantly reduced by subtracting their 1-loop perturbative expres- 
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Figure 5: Comparison of i?(xo)-data plotted vs It jT = 2(xq — uq)/T at two different 
lattice volumes for /3 = 3.90 and the smallest available values of fj,g and 



sion up to O(a 2 go), computed in [33]. Further information on the computation of 
the four-fermion operator renormalization constant, including a numerical check on 
the absence of wrong chirality mixing, is given in Appendix A. 

In the following we will quote results obtained by employing the estimates of 
Za and Z^y^y^ from the procedure M2, as well as the very precise determination 
of Zy obtained in ref. [32]. The latter is determined by comparing (matrix elements 
of) the local current to the exactly conserved one. We immediately notice that upon 
using the procedure Ml for the evaluation of Za and Z^\ VA very similar results 
for B^ Gl are obtained, with differences typically as small as 0.002 (that we include 
in the systematic error) 8 . The values of the renormalization constants entering the 
present computation of can be found in Table 4. 

3.2.2 Continuum and chiral extrapolations 

Continuum and chiral (fi£ fi u /d) extrapolations are carried out simultaneously 
exploiting the SU(2)-xPT based fit ansatz of ref. [36] which, for our computation 

8 The consistency of the continuum limit results for the unphysical quantity -B]| GI (/i|, /t£) eval- 
uated with renormalization constants from the procedure Ml or M2 was noted in sect. 2.4. 
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with n sea = fj, e , reads 



B KMW, Ms) - (Ms) 1 + b ^ s >—f2 327T 2 /g 167T 2 /o 



+ a 2 f Q 2 D B (fa), 
(3-1) 

where ^j^t (w, As) denotes the RGI lattice estimates of the "bag parameter" at 
renormalized ^-quark mass fa and strange quark mass fi s . Note that all dimensionful 
quantities are expressed in units of /q. The dimensionless fit parameters Bff^(fa), 
b(fi s ) and Db(£l s ) are functions of /t s . The first of them represents the RGI "bag 
parameter" at fi s in the SU(2)-xPT limit (fa — > 0). 

Several comments are in order here. First, as shortly mentioned above, the 
renormalized (in MS at 2 GeV) quantities Bq and fa/d and the chiral limit pion 
decay constant fo 

B Q = 2.84(11) GeV , jX u/d = 3.5(1) MeV , f = 121.0(1) MeV , (3.2) 

as well as estimates of a at /3 = 3.8, 3.9, 4.05, 

a/o|/?=3.8 = 0.0604(16) , a/ |/j= 3 .9 = 0.0514(08) , af y =4m = 0.0409(07) , 

(3-3) 

are obtained by an analysis of the data for the mass and the decay constant of the 
it (charged) pseudoscalar meson following the lines of ref. [35] and using the results 
on the quark mass renormalization constant Z M = Zp 1 from [32]. Here we just recall 
that in this analysis based on SU(2)-xPT at NLO we set the physical scale using the 
experimental value of f w and take into account finite-size effects and 0(a 2 ) artifacts. 
As we rely on NLO SU(2)-xPT formulae, the It pseudoscalar meson data are taken 
only from a fa region (see Tables 1, 2 and 3) where the impact of possible NNLO 
corrections has been checked (via the methods of ref. [35]) to be negligible. 

The renormalized strange quark mass jl s is extracted by analysing our data on 
the hi pseudoscalar meson mass M34 at (3 = 4.05, 3.9 and 3.8 with the help of the 
SU(2) xPT based ansatz [36] 



1 + C(jih) 72 
JO 



+ a 2 f%D M (fa), (3.4) 



for three reference values of fa in the strange quark mass region (from 75 to 105 
MeV). In the fit ansatz (3.4), with all dimensionful quantities expressed in units of fo, 
we have included a term, with coefficient DM(fa), to parameterize the leading 0(a 2 ) 
artifacts. Finite size effects on the raw lattice data for M| 4 , although not shown 
in eq. (3.4), were taken into account in our actual fits. These effects, which are 



9 With respect to the formulae of ref. [36] we replace A x with 47r/o and reabsorbe the effect of 
this change in a redefinition of b(fi h ). We use fo ~ 121 MeV as found in [35]. 
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checked at j3 = 3.9, afj,£ = 0.0040 to be < 1% and only marginally significant within 
errors (see Table 2), have been estimated by using resummed x~PT formulae [37] and 
subtracted out from the lattice data. The coefficients Cm (Aft) and -Dm (Ah) grow 
almost linearly with Aft in the aforementioned range, in line with the chiral behaviour 
expected at finite lattice spacing, namely M| 4 ~ (fif + /%) as (fi£ + Ah) — > [31]. 

Since at this stage the value of Au/d is already known, once the parameters 
Cm (Aft)) c(pLh) and -Dm (Aft) have been determined by the fit of the A£-dependence 
based on eq. (3.4), one immediately yields continuum limit estimates of M| 4 (A„/d, Aft) 
for the three chosen reference values of Aft- These estimates exhibit a smooth 
dependence on Aft an d can be interpolated linearly to the experimental value of 

= (495MeV) 2 , thereby providing a well controlled result for As, for which we 

find 

A s = /4 MS ' 2 GeV) = 92(5) MeV . (3.5) 

More details on this analysis will be given in a forthcoming publication [38]. 

Also our lattice estimators for B^f^ t (jli, Aft) turn out to have a weak dependence 
on Aft- It is thus straightforward to obtain by linear interpolation in Aft their values 
at the point Aft = As- By exploiting the available determinations of a(/3) in physical 
units from the analysis of the pion sector data, the interpolation of B^fl t (fif , fih) 
to the s-quark mass point can be performed separately at each value of j3 and a£i£. 
An example of such interpolations is shown in Fig. 6. In this way we obtain the 
quantity -B^^ t (A^, As) for the three lattice spacing and the aA^-values corresponding 
to the bare parameters in Tables 1, 2 and 3. This set of numbers can be fit to 
the formula (3.1) where the quantities -Bj^^As), b(fl s ) an ^ -D_b(As) are the free 
parameters to be determined. As it clearly appears from Fig. 7 the quality of the 
resulting fit is very good. 



3.2.3 Final result and error budget 

The value of B^ at the physical point is finally determined by evaluating eq. (3.1) 
in the continuum limit at the best-fit values of -B^^As) and 6(As) and setting 
V-l = fJ-u/d- I n this way we arrive at the result 

B% Gl = 0.729(25) . (3.6) 

The error quoted in parentheses in eq. (3.6) is of statistical origin. It comes from 
our fitting procedure and reflects the statistical errors on the raw data for the bare 
matrix elements and the associated renormalization constants. Its estimate has been 
performed via a standard bootstrap analysis, thereby taking properly into account 
all possible cross-correlations. As one checks in Tables 1-3, at fixed quark masses 
and lattice spacing the typical statistical error on the bare lattice estimator of Bk 
is around 1%. Inspection of Table 4 shows that relative error coming from the Bk- 
renormalization factor Zy G1 +AV / (ZaZv) is about 2.0%. The extrapolation to the 
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Figure 6: Linear interpolation of B^ Gl (fi£, (x^) to the physical strange quark mass 
point fih — > fj-s = 92 MeV at (/3 = 3.8, o/i^ = 0.0080), panel (a), (/3 = 3.9, afi£ = 
0.0040), panel (b), and (/3 = 4.05, a/j, t = 0.0030), panel (c). 



continuum limit and the physical pion mass, plus interpolation to the kaon point, 
finally leads to the error quoted in eq. (3.6). 

In order to try to have an idea of the size of the systematic error other analyses 
have been performed besides the one discussed above. 

i) First of all, the whole analysis was repeated on the same set of raw data 
but choosing a/ro extrapolated to the chiral limit [35] (with ro the Sommer scale), 
instead of afo, as the scaling variable used to build dimensionless quantities suited 
for continuum extrapolation. This change turned out to produce an increase of the 
value of B^ Gl of a few permille. 

ii) Secondly we replaced 2B^i by the squared It- meson mass (Mj™) 2 in ex- 
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Figure 7: Combined chiral and continuum fits according to the ansatz (3.1). The 
label M2 refers to the specific RI-MOM procedure adopted for the computation of 



trapolating to the physical pion mass point and we determined the kaon point by 
interpolating the quantity (Mj^) 2 to the physical value of 2M^ — M 2 . As before 
^rgi i ncreases by a f ew permille. Moreover we have tried a fit function which 
also includes higher order analytical terms added to the SU(2)— %PT formula. We 
observe a shift equal to 0.007 in the central value of B^ Gl . 

iii) We also tried to perform the extrapolation to the physical pion mass point 
using a first or second order polynomial in [/,£ in alternative to the »S[/(2)xPT-based 
ansatz (3.1). We obtained good quality fits where the central value of B^ G1 gets 
shifted by +0.021 or +0.001, respectively. Completely analogous results are obtained 
by employing a first or second order polynomial in (M^ 11 ) 2 . 

iv) Finally all these alternative analyses were repeated by using the Ml method 
for the RI-MOM renormalisation in place of M2. In all cases this produced tiny 
changes of B^ Gl (never larger than ±0.004). For instance, for the main analysis 
discussed in some detail in the present section, a decrease of B^ Gl by 0.002 was 
observed. Further information on the analyses for the estimate of the systematic 
error is given in Appendix B. 

In this way we arrive at the systematic error budget for B^ Gl we anticipated in 
the Introduction, which is detailed below. We get 

• 0.004 from the RI-MOM renormalization, evaluated by combining the uncer- 
tainty due to the spread coming from methods Ml and M2 and the estimated 
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error due the truncation to NLO of the perturbative series in the conversion 
from the RI-MOM scheme to RGI. 

• 0.009 from the uncertainty in controlling/removing cutoff effects. In this figure 
we combine the neglected 0(a 4 ) lattice artifacts which we estimate to be ~ 
1% 10 and the effect related to the choice of the scaling variable used in the 
analysis (a/o or a/ro). 

• 0.014 from the systematic uncertainty in the extrapolation to the pion mass 
point, where we add in quadrature 0.004 from varying the choice of the ex- 
trapolation variable (u^ or (M]™) 2 ), 0.011 from the difference between an 
S'C/(2)xPT-based fit and a simple polynomial ansatz (we take the average of 
the spread resulting from the polynomial fits discussed above) and a shift 
equal to 0.007 due to a possible inclusion of higher order analytical terms in 
the SU(2) X PT fit formula. . 

Summing in quadrature all these small systematic uncertainties leads to a total 
systematic error estimate of ±0.017. This number was added in quadrature to the 
statistical error in eq. (3.6), yielding the final result quoted in eq. (1.2). 

4 Conclusions 

In this paper we have computed from the Nf = 2 tm-LQCD simulation data pro- 
duced by the ETM Collaboration the value of the strong interaction parameter, Bk, 
which controls the K° — K° oscillations. To this end we have exploited the par- 
tially quenched set up proposed in ref. [22] which, at the price of introducing 0(a 2 ) 
unitarity violations, ensures O(a) improvement and absence of wrong chirality mix- 
ing effects [13]. Renormalization is carried out non-perturbatively in the RI-MOM 
scheme [25,32]. 

Using data at three lattice spacings and a number of pseudoscalar masses in 
the interval 280 MeV < mps < 550 MeV, we get in the continuum limit of Nf = 2 
QCD and at the physical value of the pion and kaon mass the value of the RGI "bag 
parameter" reported in eq. (1.2), namely B% GI = 0.729 ± 0.030. 

Our determination is quite accurate, with a ~ 4% total error (with statistical 
and total systematic errors added in quadrature), and agrees rather well with other 
existing values. In particular, from the comparison of our result for Bk with the re- 
sults from Nf = 2 + 1 dynamical simulations (see Fig. 1), it may be inferred that the 
quenching of the strange quark leads to an error which is rather small i.e. smaller, 
at present, than other systematic uncertainties affecting current Bk estimates. If 

10 The estimate of ~ 1% as the natural size of the 0(a 4 ) corrections follows from the fact that we 
observe in our B^ G1 lattice estimators 0(a 2 ) cutoff effects at a level of ~ 10%. 
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taken at face value, our result confirms the tension between the lattice determina- 
tion (which do not seem to depend significantly on the number of dynamical quark 
flavours) and the preferred value of recent phenomenological, Standard Model based 
analyses (see refs. [3,39-42]). However some theoretical concerns remain about the 
way of estimating the long distance contributions in the relation between the strong 
interaction parameter I?k and the experimentally well measured quantity ex that 
parameterizes the (indirect) CP violation in the K° — K° system. 

In order to get a substantially more accurate determination of B^ GI we would 
need to increase the statistics and work at even finer lattice spacings with the u/d- 
quark mass closer to the physical point. A further possibly important improvement 
is the use of data from more realistic simulations where dynamical strange and charm 
sea quarks are included. Actually the ETM Collaboration is already moving ahead 
in these directions [43] and work aiming at more precise lattice determinations of 
the kaon "bag parameter" is in progress. 

Acknowledgments - We wish to thank all the other members of ETMC for 
their interest in this work and a most enjoyable and fruitful collaboration. This 
work was supported by the European Union Sixth Framework Programme under 
contract MTRN-CT-2006-035482, FlaviaNet. V. G. thanks the MICINN (Spain) for 
partial support under grant FPA2008-03373 and the Generalitat Valenciana (Spain) 
for partial support under grant GVPROMETEO2009-128. F. M. acknowledges fi- 
nancial support from projects FPA2007-66665, 2009SGR502 and from the Spanish 
Consolider-Ingenio 2010 Programme CPAN (CSD2007-00042). M. P. acknowledges 
financial support by a Marie Curie European Reintegration Grant of the 7th Eu- 
ropean Community Framework Programme under contract number PERG05-GA- 
2009-249309. A.S. acknowledges financial support from the Research Promotion 
Foundation of Cyprus (Proposal Nr: ENISX/0506/17). 



Appendix A - The RI-MOM computation of Z VA+AV 

Here we discuss the non-perturbative computation of the renormalization constant 
(RC) of the operator Qvv+AA (see eq. (2.3)), which, for the reasons explained in 
sect. 2.1, we call Zva+av- It is first evaluated in the RI'-MOM scheme n , at 
some appropriate scale fio lying in the domain of applicability of (RG-improved) 
perturbation theory, following the strategy detailed in ref. [32]. Then Zy l A+AV is 

11 The prime in the label RI' is to remind the specific defining condition adopted for the quark field 
RC Z q , which, as convenient in lattice studies, is taken here to coincide with the quark propagator 
form factor Ei. For the latter we adopt the definition of eq. (33) of ref. [32]. 
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converted into the (scheme- and scale-independent) quantity Zy^^ AV according to 
the formula 



Z^aUv = («(*)))"* exp { J^ 0) da' + ^7] }zFa + av(A) > 

(A.l) 

where 

(3(a) = -|V + 0(a 3 ), ft = 11-^, 

TvawW = ?a + 0(a 2 ), 7o = 4 , (A.2) 

47T 

which is here evaluated with Nf = 2 and to NLO accuracy [44,45]. RGI quantities 
are normalized following the conventions of ref. [1]. In this appendix we quote results 
for Zy A ^ +A y but, since the conversion from RI'-MOM to RGI is a straigtforward 
step, most of our discussion will actually concern the evaluation of Zy^^dj^) 
at the three lattice spacings of interest for the present work. The error on 
associated to the NLO perturbative conversion from RI'-MOM to RGI is estimated 
to be about 0.002 and has been taken into account in the error budget discussion of 
sect. 3.2.3. This estimate is obtained assuming that the neglected NNLO and higher 
order relative corrections to the conversion factor are as large as the square of the 
NLO relative correction, which turns out to be ~ 5%. 

For the purpose of determining Zy l A+Ay , besides the gj-quark propagator 12 

VP) = a4 E e ~^M^7(0)}, (A.3) 

X 

we compute the relevant Green function is the momentum space correlator 
G vv+AA (p,p,p,p)f$ s = «"'V, e -*(*i-* 2 +*3-*4) . 

• Mx2)] b p Qvv + aa(0) Mxz)]' Mx4)] d s ) ■ (A.4) 

The where lower-case Greek (Latin) symbols denote uncontracted spin (colour) in- 



dices. In terms of the corresponding amputated Green function 13 

A vv+AA (p,p,p,p;fi val ,fi sea )f$ 5 = [S^YXi* [S g ,3(p)~%y " 

• G vv+AA (p,p,p, P r^%, [S q , 2 {p)-%^ [V(P) _1 ]™ > ( A - 5 ) 

the RI'-MOM RC is obtained by imposing in the chiral limit (fi va i = Usea = 0) the 
renormalization condition 

[ Z VA+AV Z q 2 D n](P] Hval , Usea) = 1 , (A.6) 



12 The equations below are written in the so-called physical quark basis with the lattice valence 
quark action given by eq. (2.2). 

13 We explicitly display in the l.h.s. of the equation the dependence on quark mass parameters 
that is implicit in the quantities appearing in the r.h.s. 
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where Z q denotes the (flavour independent) quark field RC and the VV + AA- vertex 
D^i is given by 

Dn(p; fi V al,Vsea) = Aw+Aa(P, P, P, V, Vval , ^sea)afa' P'^IpJ^ ( A - 7 ) 

where V^^ A is a suitably normalized spin-projector reflecting the VV + AA Dirac 
structure of the vertex. Details on this projector and the notation can be found 
in [15]. 

As usual in lattice RI-MOM computations, the relevant Green functions (A. 3)- 
(A.4) are evaluated in the Landau gauge for a sequence of sea {fJ, S ea} and valence 
{^vai} quark mass parameters at each of the three lattice spacings we consider 
here. The bare parameters and the statistics of this computation are detailed in 
the Table 2 of ref. [32]. Also the lattice momenta p u , with pi,2,3 = (27r/-L)ni 5 2,3, 
P4 = (27r/L)(n4 + 1/2), we consider here are specified in eqs. (20)-(21) of ref. [32]. 
To minimise the contributions of Lorentz non-invariant discretization effects we have 
considered in our analysis only the momenta p satisfying the constraint 

P$ < 0-28 (J2 fu? , aPu = sm(aPu) ■ (A.8) 

P v 

In the following we shall often use the quantity p 2 = J2vft1- 



Sketch of procedure for extracting Zy A+AV 

Here we summarise the rather standard procedure we follow to extract Zy A+AV , 
deferring the discussion of further details to subsequent subsections. First of all, for 
each j3 and choice of the scale p 2 , the condition (A. 6) is enforced at all the values of 
fi sea and fi va i given in Table 2 of ref. [32]. By doing so one obtains lattice approx- 
imants at non-zero quark mass(es) of the desired RC, Zy l A+AV (p 2 ; a 2 p 2 ; n va u ^sea) 
which are then extrapolated to fi va i = fi sea = 14 . 

Improved and controlled estimates of Zy A+AV (p 2 ; a 2 p 2 ; 0,0) are obtained by 
removing the perturbatively leading cutoff effects. Using NLO continuum QCD 
evolution [44,45], the first argument is brought to a reference scale value, /j,q. The 
reliability of this step rests on the usual assumption that the scales p 2 and are 
large enough to make NLO-PT accurate (within statistical errors). The remaining 
a 2 p 2 dependence in Zy l A+AV ($ \ a 2 p 2 ; 0, 0) is to be regarded as a mere lattice artifact 
which will be taken care of by employing either the Ml or M2 method (introduced 
in ref. [32]), as briefly described below. 

In order to reduce the statistical error, the lattice approximants of the Qvv+AA 
RC have been averaged over equivalent pattern of Wilson parameters (ri, r2, r%, r±) 

14 To facilitate the discussion, in the notation for these quantities we shall distinguish, as it is 
customary (see e.g. [32]), the scale-dependence described by the RG eq's of continuum QCD (first 
argument) from the one due to 0(a 2 p 2 ) cutoff effects (second argument). 
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p 




■^va+av (M2) 


Z A {Ml) 


Z A (M2) 


Z V 


3.8 


0.605(18) 


0.620(13) 


0.746(11) 


0.727(07) 


0.5816(02) 


3.9 


0.623(11) 


0.630(08) 


0.746(06) 


0.730(03) 


0.6103(03) 


4.05 


0.691(11) 


0.698(08) 


0.772(06) 


0.758(04) 


0.6451(03) 



Table 4: Results for Zy^ j4 y(Ml,M2) from the analysis discussed in this appendix, 
together with the (best) final estimates of Za{M1,M2) and Zy from ref. [32], for the 
three /3-values of interest here. 

of valence quarks (namely (1, 1, 1, —1) and (—1,-1,-1,1)) as well as over different 
lattice momenta corresponding to the same value of p 2 . We have checked that 
performing these averages before or after taking the chiral limit leads to perfectly 
consistent (actually almost identical) results. The outcome of the whole analysis is 
conveniently expressed in terms of Zy^ Av (Ml) and Zy A iI +AV (M2). These are the 
RGI quantities whose values are quoted in Table 4 together with the values of the 
other RC's relevant for the B^ Gl computation. 

It is interesting to compare our results with those obtained from a 1-loop pertur- 
bative lattice computation. For instance in the case /3 = 3.9 one gets Zy'^^y" 01 = 
0.422, if the perturbative gauge coupling is set to g 2 ./ (P) ((P) is the plaquette ex- 
pectation value), or Z v ~ A l( ^ > A y' Gl = 0.675, if the perturbative gauge coupling is set to 
the "boosted" value corresponding to the prescription of ref. [46]. These numbers 
should be compared with the values in Table 4. 



Valence chiral limit 

For fixed values of f3, a 2 p 2 and a/j, sea we fit D\\ (see eq. (A. 7)) to the ansatz 

Du(p; Hvah Usea) = A(p 2 -{l sca ) + B(f- /J, sea ) /J, va l + C(p 2 ] /J, sea ) (/J-val)' 1 , (A.9) 

the form of which we would like now to justify. 

In general the RI-MOM approach relies on the fact that, at very large values of 
p 2 (S> Aq CD ), the relevant Green functions (and the associated vertices) are poly- 
nomial in all the quark mass parameters, because non-perturbative contributions 
which can violate this property vanish in this limit. However, at finite values of p 2 
special care must be taken in studying the quark mass behaviour of Green functions 
that admit one Goldstone boson (GB) intermediate states in their spectral decom- 
position. Indeed, such non-perturbative contributions to Green functions, though 
suppressed by a factor 1/p 2 for each GB pole (see below), are divergent in the chiral 
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Fit D n (.p;fi v;iJ .(i s|;a )=A+Bn val +C/(i vji| 
GB-pole subtracted D 



Fit D n foji v .,|r^=A+Bu, +<Vu, 
GB-pole subtracted D 




(a) 



(b) 



Figure 8: GB-pole subtraction and valence chiral limit of D\\ (see eq. (A. 9)), here 
plotted vs. aii va t, for f} = 3.9, an sea = 0.0040, a 2 p 2 ~ 1.565 (left panel) and /3 = 4.05, 
a^sea = 0.0030, a 2 p 2 ~ 1.568 (right panel). 



limit and must hence be disentangled and removed. As explained in Appendix B of 
ref. [47], this is precisely the case for the correlator Gyv+AA (see eq. (A. 4)) and the 
associated vertex D\\. 

By using the LSZ reduction formulae one finds that the GB pole contributions 
relevant for (the spectral decomposition of) Du(p; fj, va i, fj, sea ) are of the form 15 

W ^^ (P M \Qvv + AAqi(p)q2(-p)\n) P 2 , 

(h) { M?2y p 2 (%3fr)g4(-p) Qvv + aa\P 12 )p 2 , 
A 2 A 2 

, s A QCD /n 34i^ | „12\ QCD 

H (M 34)2 p 2( P \Qw+Aa\P ) {Ml2) 2 p 2 > 

where, for simplicity, all the coefficients that are regular in the GB squared mass 
and weakly (at most logarithmically) depending on p 2 have been set to unity. Fur- 
thermore immaterial (for this discussion) lattice artifact corrections are neglected. 
As discussed in sect. 2, owing to the choice (2.4) of the valence quark Wilson param- 
eters, the axial current 947^75(73 is exactly conserved on the lattice (only broken by 
soft mass terms), while the conservation of (?27/i75<Zi is spoiled by lattice artifacts. 



15 Here qi(p) denotes the four-dimensional Fourier transform of qi(x) with momentum p. As in 
sect. 2 M 12 (M 34 ) is the mass of the lattice GB state |P 12 } (|P 34 )). 
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From the form of Qw+AA it follows in particular that the matrix elements 

(P 34 \Q VV+AAqi (p)q 2 (-p)\n) , (P 34 \Q VV+AA \P 12 ) 

vanish as (M 34 ) 2 ~ (/X3 + ^4) in the limit ^4 — > 0. This property can be checked by 
using for these matrix elements the soft pion theorem associated to the conservation 
of the lattice current <?4 7^7593 in order to reduce the GB-particle P 34 . Hence in the 
spectral decomposition of D\\ no terms with two GB-poles occur. Moreover, the 
terms (ii) and (iii) above, which both contain one GB-pole, are strongly suppressed 
by factors of the kind a 2 /p 2 and (1/p 2 ) 2 , respectively. 

In line with these arguments a very smooth dependence of D\\ on /_t va i (or 
equivalently on M 2 2 ) is observed for all the momenta of interest in our analysis, 
see Fig. 8 for two typical examples. By fitting the D\\ data to the ansatz (A. 9) 
the valence chiral limit is thus safely determined. A practically identical result is 
obtained by employing a fit ansatz analogous to (A. 9), with [i va i substituted by M 2 2 . 
Combining the valence chiral limit lattice estimator of Dn with the corresponding 
approximant of Z q (the valence chiral extrapolation of which poses no problems [32]) 
leads to reliable estimates of the intermediate quantities Zy A+AV (p 2 ; a 2 p 2 ; 0, a/j, sea ). 




Figure 9: The quantity Zy l A+AV (p 2 ; a 2 p 2 ;0, afi sea ), taken at the valence chiral limit, 
as a function of a 2 fi 2 ea , for a typical lattice momentum choice (see inset) giving 
a 2 p 2 ~ 1.06, for the three j3 values considered in this paper. 
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Sea chiral limit 



This limit is taken at fixed f3 and a 2 p 2 by fitting the Zy l A+Ay {n^] a 2 p 2 \ 0, a/x sea )-data 
to a polynomial of first order in a 2 fj, 2 ea . Within our statistical errors (typically at 
the level of 1 to 2%) on the data at fixed a(j, sea , the dependence on a// sea is hardly 
significant, as one can see e.g. from Fig. 9. We checked that repeating the whole 
analysis with Zy A+AV (fiQ; a 2 p 2 ; 0, a/j, sea ) fitted instead to a constant in afi sea leads 
to compatible results for the desired RC, but with smaller errors and often (but 
not always, in particular at j3 = 3.9) acceptable x 2_va lues. Based on these findings 
we conservatively decided to perform the sea chiral extrapolation via a linear fit in 
a 2 fi 2 ea . In this way we get the RC-estimators Zy A+A y(/j,Q] a 2 p 2 ; 0, 0). 




Figure 10: The effect of subtracting from Z$ I A+Av (p 2 ;a 2 p 2 ;0,0) at /3 = 4.05 (blue 
dots) the 0(a 2 g 2 ) correction (A. 10), setting either g 2 = (red squares) or g 2 = g 2 
(green diamonds). 
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Removal of 0(a g ) cutoff effects 

Improved chiral limit RC-estimators, Zy l A l { _ 1 A y r (p 2 ;a 2 p 2 ), are obtained by removing 
from Zy l A+AV {p 2 ] a 2 p 2 ; 0, 0) discretization errors up to 0(a 2 <7 2 ), viz. 

<;7v P V; «V) = Z^ +AV (p 2 ; a 2 p 2 ; 0, 0) 



9 2 2 



p 2 (4 i ; 1 +^: 1 iog(aV))+ c ^ l 



»(2) 



,(3) 



.9 



12vr 2 



2 a 2 



~2l (!) , ( 2 ) t / 2~2\\ , (3)^pPp 
P {Cq ' +C q ' log(aV)) + 4 ^ 2 



(A.10) 



The coefficients c, 



C?) „■ 

9 ) J 



1,2,3 can be found in eq. (34) of ref. [32], while the 



coefficients , j = 1,2,3, namely 



2.64223 , 



J2) 
-Da 



-19/18, 



" ( n 3) 



-2.79899, 



(A.11) 



have been recently computed by evaluating the vertex Dn{p) (eq. (A. 7)) in lattice 
perturbation theory, to 1-loop and including cutoff effects up to the second order in 
a [33]. Here we just recall that in lattice perturbation theory the vertex Dn(p) can 
be expressed in the form 



£>n(p) pert = 1 + 



g 



16vr 2 



C + b ( ul log(a 2 p 2 )) 



Wir 2 



p 2 (cg; i +c^ 1 log(aV)) + c^ 1 



.(2) 



,(3) 



+ 0(aV,5 4 ) (A.12) 



(1 2) 

and refer to [33] for the values of the (here irrelevant) coefficients b D ' and all details. 
The form of eq. (A. 10) follows directly from eqs. (A. 6) and (A.12) above, as well as 
from eq. (32) of [32], where the perturbative expression of = Z q (p), including 

the 0(a 2 g 2 ) corrections, is given. 

In the numerical evaluation of the perturbative correction in eq. (A. 10) g 2 
was taken as the simple boosted coupling g 2 = g^j (P) , where the average pla- 
quette (P) is computed non-perturbatively. The values of (P) we employed here 
are [0.5689,0.5825,0.6014] for /3 = [3.8,3,9,4.05], respectively. The (nice) effect of 
the (A. 10) correction in removing the unwanted a 2 p 2 dependence is illustrated, in the 
case (3 = 4.05, in Fig. 10. In the figure the uncorrected values of Zy A+AV (p 2 ; a 2 p 2 ;0, 0) 

are compared with the values of Zy A ~^ A y T (p 2 ; a 2 p 2 ) obtained by setting either g 2 

9l ^ 



— 6//3 or (as we chose in the end) g 2 = g 2 . 
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Absence of wrong chirality mixings 

According to the analysis of ref. [22], at maximal twist wrong chirality mixings can 
affect the renormalization of Qvv+AA only at order a? (or higher). The relation 
between the renormalized lattice and bare operators has then the form 

5 

QvV+AA = Z VA+AV [Qw+AA + ^ &ljQj + •••], (A. 13) 

J'=2 

where (using the operator basis of ref. [15]) we have 

Q2 = 2 {[<Zi7m<?2] [(737/^4] - [917^7592] [<737m7594] + (92 <-> 94)} , 

Q 3 = 2{[q 1 q 2 ][q- i qA\ - [917592] [937594] + (92 94)} , 

Q 4 = 2{[gig 2 ]fe94] + [9175^2] [937594] + (92 <-> 94)} , 

Q 5 = 2 {[910-^92] [93^94] + (92 ^94)} • ( A - 14 ) 

In eq. (A. 13) the mixing coefficients A±j are 0(a 2 ) quantities and . . . stand for 
higher order lattice corrections. In Fig. 11 we plot Ay, j = 2,3,4,5 as function 
of a 2 p 2 for (3 = 3.8, /3 = 3.9 and j3 = 4.05. It is clearly seen that in all cases the 
mixing coefficients are zero within errors and that this is systematically more so as 
/3 increases. 



Final estimates from Ml and M2 methods 

The first step is to convert Zy]^ 1 ^ (p 2 -,a 2 p 2 ) to Z v 1 A ^^y C {^; a 2 p 2 ) by using the 
known formula for the NLO running [44,45]. This step is necessary in order to 
disentangle the 0(a 2 p 2 ) cutoff effects from the genuine continuum p 2 dependence, 
but the actual value of fio has no impact on the RGI result for the RC. As customary, 
we take ^0 = V a (/^) f° r eacn value of (3, namely l/a(3.8) = 2.0 GeV, l/a(3.9) = 
2.3 GeV and l/o(4.05) = 2.9 GeV. 

The method Ml consists in fitting Z^" 1 ™^^; a 2 p 2 ) to the ansatz 

Z^A+Av^h a V) = ^VA+Avi^o) + ^VA+Ava 2 p 2 (A. 15) 

in the large momentum region, 1.0 < a 2 p 2 < 2.2. As expected, the slope Xva+av 
exhibits a very mild gauge coupling dependence which (unlike the case of quark 
bilinear RC's) is not significant within our current statistical errors. We thus treat 
X-va+av as a /3-independent quantity. According to the ansatz (A. 15) a simultaneous 
linear extrapolation to a 2 p 2 = at the three f3 values of interest was performed. The 
extrapolated values, Zy A . av {^q), are finally used to evaluate via eq. (A.l) to NLO 
accuracy the quantities 

Z va+av( M1 ) quoted in Table 4. 
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Figure 11: The behaviour of the mixing coefficients Ay, j = 2, 3, 4, 5 as function of 
a 2p2 £ Qr p _ 3 p _ an( j p = 4 05) reS p ec tively. 



The numbers for Zy 1 ^^^ (hq] a 2 p 2 ) at the three /3's and the lines of simulta- 
neous best linear fit in a 2 p are shown in Fig. 12. We recall that in this figure, as 
in all other figures presented in this Appendix, only data points corresponding to 
the momenta p satisfying the constraint (A. 8) appear. As we said above, data that 
refer to momenta p violating the condition (A. 8) are never used in our RC analysis. 

When the method M2 is used, we separately average at each /3 the values 
of Zy l j^^y T a 2 p 2 ) over a narrow interval of momenta where data are linear in 
a 2 p 2 . In physical units we take the same momentum interval for all /3's, i.e. p 2 6 
[8.0,9.5] GeV 2 . 
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Figure 12: Z^/^ r (a(/3) 2 ;a 2 p 2 ) as a function of a 2 p 2 for the three f3 values con- 
sidered in this paper. The three straight lines represent the simultaneous linear fit 
to the Zy l ' A ~^\a(l3)- 2 ; a 2 p 2 ) data in the interval 1.0 < a 2 p 2 < 2.2 at the three /3's. 



Appendix B — Chiral and continuum extrapolations 

Here we give some details and discuss a few technical aspects of the auxiliary analyses 
that were performed, as discussed in sect. 3.2.3, to estimate the systematic error on 
our Bft G1 computation. 

In order to evaluate the systematic errors affecting the extrapolation in the 
u/d-quark mass, which was performed together with that of the continuum limit, 
we also tried an analysis where the role of the renormalized quark masses, p,£ and jl^, 
is played by (M^™) 2 and (M^ 1 ) 2 i.e. the masses of the pseudoscalar mesons made 
out of two mass-degenerate quarks (regularized with opposite Wilson parameters, as 
implied by the label "tm"). The fit ansatz (3.1) was then replaced by the following 
one 

Bl%{M\?,M%{)= (B.l) 



i +h '{M^y—^-"—^\og 



f 2 32vr 2 / 2 & 16vr 2 / 2 



+ a 2 f 2 D' B (M^) 



where B'^ l (M^), b'(M^) and D' B (M}%) are the three fit parameters. For all 
values of /3 and M^ m , data for B^ &t {Mjf,M^) have been first interpolated in 
(M^™) 2 16 in order to obtain the RGI bag parameter at three reference values of 



16 The interpolation was based on global fits linear in (A/£™) 2 . The fit quality was good in all cases 
and the statistical error on the interpolation results has been estimated by a standard bootstrap 
procedure. 
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Figure 13: Combined chiral and continuum fits according to the ansatz (B.l) yielding 
B^ GI at fi h corresponding to Mffi/f = 5.50. The label M2 refers to the RI-MOM 
procedure adopted for the computation of ZyA+AV anc ^ 



Mfrfi, namely 5.10/o, 5.50/o, 5.90/o- Then three continuum and chiral fits for B^i^ 
at these three reference values of Mj^ were performed, based on the ansatz (B.l). 
The fits turned out to be of good quality and yielded continuum limit results for 
B% GI (M n ,M hh ) at M n = 135 MeV and M hh = (5.10, 5.50, 5.90) / , respectively. For 
illustration we report in Fig. 13 the results for the case of Mj^ = 5.50/o- The 
extrapolated value, Bg (M^, 5.50/o), is represented by an open circle in the figure. 
As one can see from Fig. 14, the dependence of these extrapolated results on Mf ih is 
very mild and hardly significant compared to the 4%-level of our statistical errors. 
In this situation the interpolation of B^ G1 (M 7r , Mhh) in (-^™) 2 to the physical 
kaon mass (which we identified as the point where (M^ 1 ) 2 = 2M 2 K - Ml) poses no 
particular problem (see Fig. 14), and yields BW* 1 = 0.732(24). Notice that owing 
to the weak dependence of B^ (M^, Mhh) on (Mhh) 2 ' the approximation induced 
by the use of the LO SU(3)-xPT formula to relate {Miff to the physical value of 
the 2M^ — M 2 mass difference turns out to have a completely negligible impact on 
the final result. 

With the purpose of estimating the systematic uncertainty involved into the 
extrapolation of our Bx data to the pion (or -u/d-quark) physical point, we also 
studied the effect of performing the chiral fit assuming a first or second order poly- 
nomial dependence on fj,g (or (M^™) 2 ). The analysis of the fah and quark mass 
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Figure 14: B^ Gl at the "physical" u/d quark mass as a function of (M^™//o) 2 . 
Squares represent the values computed as explained in the text at M^h = 
(5.10, 5.50, 5.90)/ . The interpolated value at the physical point (M^) 2 = 2M\ - 
M 2 is indicated with a full dot. The label M2 has the same meaning as in Fig. 13. 



dependence, described in sect. 3.2.2, was repeated with the chiral and continuum fit 
ansatz (3.1) replaced by 

« t (A^ s ) = ^ G x I (/i s )+C 1 (/i s )^+C 2 (A s )f^') +a 2 / 2 C L (A s ), (B.2) 

Jo \ Jo ) 

while always taking jj, s = 92(5) MeV. In the case where only a linear dependence 
on fj,£ is assumed, we set C2(/u s ) = 0. Besides B^ G ^(fi s ), the other fit parameters 
in the formula (B.2) are Cz,(/t s ), Ci(/t s ) and (only for the case of a second order 
polynomial) C^As)- 

In the case of an assumed quadratic dependence of BW^L data on fig, in the con- 
tinuum limit and at the physical u/d-qu&vk mass point one finds B^ Gl = 0.730(51). 
The central value of this determination is perfectly in line with all our previous 
evaluations, only the associated statistical error is somewhat larger because here we 
have been dealing with a four (instead of a three) parameter fit. The good quality 
of the fit is well illustrated by Fig. 15. A good fit yielding B^ Gl = 0.750(26) is 
also obtained if just a simple linear dependence of the data on /t£ is assumed. Not 
unexpectedly (given the distribution of data points) the central value of B^ Gl in this 
case turns out to be slightly larger than in chiral fits that are based on SU(2) xPT 
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or an ansatz allowing for some curvature in the (i£ behaviour. This effect, which is 
hardly larger than one (statistical) standard deviation, is taken into account in our 
systematic error budget (see sect. 3.2.3). 
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Figure 15: Combined chiral and continuum fits B^ according to the ansatz (B.2). 
The label M2 has the same meaning as in Fig. 13. For comparison the continuum 
limit curve, corresponding to the ansatz (3.1) based on SU(2) is a l so shown. 



For completeness we also show in Fig. 16, the chiral and continuum extrap- 
olation of our Bft Gl data when the method Ml instead of M2 (see ref. [32] and 
Appendix A) is employed in the RI-MOM computation of Za and Zy^+AV By an 
analysis otherwise identical to that discussed in sect. 3.2.2 we obtain in the contin- 
uum limit and at the physical pion mass point B^jP* = 0.727(30). As expected from 
the RC values reported in Table 4, in the present case the cutoff effects at finite 
lattice spacing are somewhat smaller but the statistical errors a bit larger compared 
to the case of Fig. 7. 

As stated in sect. 3.2.3, the effect of replacing with a/ro the scaling variable 
a/o to build dimensionless quantities in intermediate steps of our analyses and mea- 
sure the distance from the continuum limit, is as small as 0.004-0.005. The curves 
illustrating our interpolation to the s-quark mass point and extrapolation to the 
u/d-quark point and continuum limit in the case where a/ro is used as scaling vari- 
able look almost identical (within graphical resolution) to those we presented above 
and are hence not shown. 

Other analyses were also performed corresponding to different variants of our 
reference analysis, discussed in sect. 3.2.3. In all cases the values of B^ Gl at the 
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Figure 16: Combined chiral and continuum fits of BjjP according to the ansatz (3.1), 
for the case where the Ml method is adopted for the RI-MOM computation of both 
-^FA+AV an d Za ~ as implied by the label Ml in the ordinate. 



physical point were found to be very close to the value quoted in eq. (3.6), with a 
spread that is well covered by our final estimate, ±0.017, of the systematic error. 
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